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The Boltzmann equation for inelastic Maxwell models is considered to determine the velocity 
moments through fourth degree in the simple shear flow state. First, the rheological properties 
(which are related to the second-degree velocity moments) are exactly evaluated in terms of the 
coefficient of restitution a and the (reduced) shear rate a*. For a given value of cv, the above 
transport properties decrease with increasing shear rate. Moreover, as expected, the third-degree 
and the asymmetric fourth-degree moments vanish in the long time limit when they are scaled with 
' the thermal speed. On the other hand, as in the case of elastic collisions, our results show that, 

for a given value of a, the scaled symmetric fourth-degree moments diverge in time for shear rates 
larger than a certain critical value a* (a) which decreases with increasing dissipation. The explicit 
shear-rate dependence of the fourth-degree moments below this critical value is also obtained. 
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I. INTRODUCTION 

o 

On of the most widely studied inhomogcncous states in granular gases is the so-called simple or uniform shear flow 
(USF) [l], This state is characterized by a constant density n, a uniform granular temperature T, and a linear 
velocity profile u x = ay, where a is the constant shear rate. The temperature changes in time due to two competing 
effects: the viscous heating and the inelastic collisional cooling. Depending on the initial condition, one of the effects 
, prevails over the other one so that the temperature either increases or decreases in time, until a steady state is reached 
for sufficiently long times. After a short kinetic regime, the time evolution and the steady state of the system admits 
a non-Newtonian hydrodynamic description [3, |4j characterized by shear-rate dependent viscosity and normal stress 
, O, ' differences. 

The prototypical model of granular gases consists of inelastic hard spheres (fHS) with a constant coefficient of 
£S) . normal restitution a < 1. fn the low-density limit, all the relevant information on the system is provided by the 
' one-particle velocity distribution function /(r, v;t), which obeys the Boltzmann equation [5[. However, it is generally 
not possible to get exact analytical results from the Boltzmann equation for fHS, especially in far from equilibrium 
situations such as the USF. Consequently, most of the analytical results reported in the literature have been derived 
by using approximations and/or kinetic models @, 0, S, [3, El El El El EI EE El E3- 

The lack of exact analytical results can be overcome in some situations by considering the so-called inelastic Maxwell 
models (IMM), where the collision rate is independent of the relative velocity of the two colliding particles. These 
models have received a lot of attention in the last few years since they allow for the derivation of a number of exact 

results [n, ei m m m m m m m m m shhiiiiii m m, m m, m no, no, m m, in . Therefore, 

the influence of inelasticity on the dynamic properties can be analyzed without introducing additional, and sometimes 
uncontrolled, approximations. In addition, it is interesting to remark that recent experiments [4o| for magnetic grains 
k> , with dipolar interactions turn out to be well described by IMM. 

In the context of the USF, the rheological properties, which are related with the second-degree velocity moments, 
have been obtained exactly in the steady state from the Boltzmann equation for IMM [2l|, [37[. However, even though 
these properties are physically important, they provide a partial piece of information about the velocity distribution 
function /, especially in the high-velocity region, where higher degree velocity moments play a prominent role. By 
symmetry reasons, the third-degree moments arc expected to vanish in the USF. Therefore, the first non-trivial 
moments beyond the ones associated with the rheological properties are the fourth-degree moments. Their knowledge 
provides relevant information about the combined effect of shearing and inelasticity on the velocity distribution. 

The effort of going from second-degree to fourth-degree moments in the USF problem can be justified by a number 
of reasons. For instance, their knowledge is needed to evaluate transport properties in situations slightly perturbed 
from the USF state 46], which allows one to perform a linear stability analysis around that state. Another interesting 
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issue is to explore whether or not the divergence of the fourth-degree moments for elastic Maxwell molecules beyond 
a certain critical shear rate [47l. l48l. |49| is also present in the inelastic case and, if so, to what extent. 

The main aim of this paper is to determine the fourth-degree moments of IMM subject to USF. This can be carried 
out thanks to recent derivations by the authors of the fourth-degree collisional moments for IMM [44| . Those moments 
are proportional to an effective collision frequency which in principle can be freely chosen. Here we will consider 
two classes of IMM: (a) a collision frequency vq independent of temperature (Model A) and (b) a collision frequency 
vq{T) monotonically increasing with temperature (Model B). While Model A is closer to the original model of Maxwell 
molecules for elastic gases [43,l5fJ], Model B, with vq(T) oc T 1 / 2 , is closer to IHS. The possibility of having a general 
function v(T) is akin to the class of inelastic repulsive models recently introduced by Ernst and co-workers [ill |42|. 
As will be shown below, Model A and B yield the same results in the steady state. In particular, the reduced shear 
rate a* = a/v in the steady state is a "universal" well-defined function a* (a) of the coefficient of restitution a. The 
main feature of Model A is that a* does not change in time and so a steady state does not exist, except for the 
specific value a* = a* (a). However, a non-Newtonian hydrodynamic regime (in which a* and a are independent and 
arbitrary parameters) is reached for asymptotically long times. This allows to study analytically the combined effect 
of both control parameters on the (scaled) velocity moments for Model A. 

The plan of the paper is as follows. In Section [Til the Boltzmann equation for IMM is introduced and the ex- 
plicit expressions for the collisional moments through fourth-degree are given. Section [Ml deals with the rheological 
properties (a nonlinear shear viscosity rj* and a viscometric function 'J) of the USF state, which are related to the 
second-degree velocity moments (pressure tensor) . While Model A lends itself to obtain the exact forms of rf and \t 
as nonlinear functions of a* and a, that is not the case for Model B since those rheological quantities require to be 
solved numerically, except in the steady state. For elastic collisions (a = 1), our expressions of 77* and "J obtained 
for Model A reduce to the results derived long time ago by Ikenberry and Truesdell [51( for Maxwell molecules. The 
third- and fourth-degree moments for Model A are analyzed in Section [TV] As expected, the results show that, when 
the third- and asymmetric fourth-degree moments are conveniently scaled with the thermal speed, they vanish in the 
long time limit. This is not the case for the symmetric fourth-degree moments. In a way similar to the case of elastic 
Maxwell molecules [U EH, El|]> we find that, for a given value of the coefficient of restitution, those moments diverge 
in time for shear rates larger than a certain critical value a* (a), which decreases as a decreases. Below this critical 
value, the (scaled) fourth-degree moments have well-defined values in the long time limit. The paper is closed in 
Section [Vl with some concluding remarks. 



II. THE BOLTZMANN EQUATION FOR IMM. COLLISIONAL MOMENTS 

In the absence of external forces, the inelastic Boltzmann equation for IMM reads 

(5 t +v-V)/(r,v;t) = J[v|/,/], (1) 
where the Boltzmann collision operator J[v|/, /] is given by [36[ 

J[vl|/ ' /] = d^/ dV2 / d? [""Wi)^) - /( v i)/( v 2)] • (2) 

Here, 

n = [ dv/(v) (3) 



is the number density, v is the collision frequency (assumed to be independent of a), Qd = 27r d / 2 /r(d/2) is the total 
solid angle in d dimensions, and a < 1 refers to the constant coefficient of restitution. In addition, the primes on the 
velocities denote the initial values {v'^vJ,} that lead to {vi,V2} following a binary collision: 

vj =v 1 -i(l + a- 1 )(?-g)5 ; , V 2 = v 2 + ^ (1 + a- 1 ) {B ■ g)a , (4) 

where g = vi — V2 is the relative velocity of the colliding pair and <x is a unit vector directed along the centers of the 
two colliding particles. Apart from n, the basic moments of / are the flow velocity 



u=-/dvv/(v) (5) 
n J 
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and the granular temperature 

T=£|dv^ 2 /(v), (6) 

where V = v u(r) is the peculiar velocity. The momentum and energy fluxes are characterized by the pressure 
tensor 



Pii = m J dvV^/(v) (7) 

and the heat flux 

q=||dvF 2 V/(v). (8) 

Finally, the rate of energy dissipated due to collisions defines the cooling rate £ as 

m 



dnT 



JdvV 2 J[v\f,f}. (9) 



The main advantage of the Boltzmann equation for Maxwell models (both elastic and inelastic) is that the (colli- 
sional) moments of J can be exactly evaluated in terms of the moments of /, without the explicit knowledge of the 
latter 50]. This property has been recently exploited [44[ to obtain the detailed expressions for all the third- and 
fourth-degree collisional moments as functions of a in d dimensions. In order to get the collisional moments, it is 
convenient to introduce the Ikenberry polynomials (50j i2r|iii 2 ...i s (V) of degree 2r + s. The Ikenberry polynomials of 
degree smaller than or equal to four are 

r | (v) = i, Y 0]i (v) = Vi, (io) 

^2|o(V) = F 2 , Y m (V) = V i V j -^V 2 S ij , (11) 
Y 2 \i(V) = V 2 V U Y om {V) = ViVjV k - j^V 2 {ViS jk + V 3 5 lk + V k 8 l0 ) , (12) 



Y M0 (V) = V 4 , F 2 |y(V)=y 2 (V i V j --V 2 S ij ), (13) 



1. 

d 

Y ome (V) = W 3 V k V e - j^V 2 {ViVjSu + W k S jt + ViVeS jk + V V k 8 lt + VjV e S ik + V k VA,) 

+ (d + 2)(d + 4) ^ ^ ij6ki + 5ik5jl + SieSjk ^ 
= ViVjVkVi - ^ [Y 2lij (V)5ki + Y 2{ik (V)5 ji + Y 2{u (V)5 jk + Y 2{jk (V)5 u 

+Y 2{ji {\)S ik + F 2 | fe ,(V)%] - d ^ +i) V A {Sij&u + 8 lk 5 3 t + 5 it 5 jk ) . (14) 

The corresponding velocity moments M2 r \i 1 i 2 ...i s and collisional moments J 2 r\i 1 i 2 ...i s arc defined, respectively, as 

M 2r \ lll2 .., s = I dvy 2r | ili2 ... is (V)/(v), (15) 



J2r\i\%2 ■ --is 



J dvF 2r | ili2 ... is (V)J[v|/,/]. (16) 



In particular, Af | = n, M ^ = (by definition of the peculiar velocity), M 2 \o — pd/m, where p = nT is the 
hydrostatic pressure, M ^j = (Pij — pSij) / 'm, and M 2 u — 2qi/m. Moreover, conservation of mass and momentum 
implies J |o = and Jg|i = 0, respectively, while J 2 \q = —C,M 2 \q- 
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The explicit expressions for the collisional moments J2r|i 1 i 2 ...i s for 2r + s < 4 are [44 

■h\0 — — C^2|0i Jo\ij — —Vo\2M \ij, 
Jl\i — — v 2\\Ml\i, Jo\ijk — — ^0\3^0\ijk} 



4|0 



v 4 \o M 4\o + Am 1 M% 10 - X 2 n L M \ ij M \ ji , 



(17) 
(18) 
(19) 



-i^| 2 M 2 |y + A 3 n 1 Ma| M |y - A 4 n 1 f M | ifc M |]y - ~M | W M |«Ai) 



(20) 



J(s\ijM — — v a\iMQ\ijM + A 5 n 



M \ijM \ M + M \ ik M y e + M \ u M \ jk - {M \ ip M \ P jdu 



+MQ\ ip M Q \ pk 8jt + M \ ip M \ p i8.j k + M \j p M \ p ),6u + M ^ p M ^ pe 5 lk + M \ kp M \ pi 8ij) 
2 



+ 



(d+2)(d + 4) 



-^0|pg-^0|qp (8ijS k £ + Si k Sj£ + SuSjk) 



(21) 



In Eqs. (|19p - (|21[) . the usual summation convention over repeated indices is assumed. The cooling rate £ and the 
effective collision frequencies f 2r i s are given by the expressions 



, d + 2 , 2 . 
C=^(l-a 2 )-o, 



(22) 



^0|2 = C + 



(l + a) s 



(23) 



3 A (l + a) 2 (d-l) 
^ = 2 C+ " 47 ^' 



^0|3 = 2^012, 



(24) 
(25) 



(1 + a) 2 (Ad - 7 + 6a - 3a 2 ) 
^4|o = 2C H — vq, 



1U 



(26) 



(1 + a) 2 [3d 2 + 7d - 14 + 3a(d + 4) - 6a 2 l 
" 2 I 2 = 2C + M{d + A) «*" 



(27) 



(1 + a) 2 [d 3 + 9d 2 + I7d - 9 + 3a(d + 4) - 3a 2 ] 
1/014 " 2C + 2d(d + A)(d + 6) "°' 



(28) 



where we have called v$ = 2v/{d + 2). According to Eqs. (jTTJ) and (f2"3")l , i/q represents the effective collision frequency 
associated with the shear viscosity in the elastic limit [52|]. Finally, the cross coefficients in Eqs. (fl"9|) - (|2"Tj) are 



(1 + af{d + 2) (Ad - 1 - 6a + 3a 2 ) 
Al = 1M 2 ^ 



(29) 



(1 + a) 2 (l + 6a-3a 2 ) 
= 8d 



(30) 
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\d 2 + 5d - 2 - 3a{d + 4) + 6a 2 l . x 
-vo, (31) 



8d 2 

(l + af [2-d + 3a(d + A)-6a 2 ] 

A4 " M(dTI) "°> (32) 

_{l + a) 2 [d 2 + 7d+9-3a{d + A) + 3a 2 ] 
h = 2d(d + 4)(d + 6) V °- (33) 

The above results hold independently of the specific form of the collision frequency isq. On physical grounds, 
i>q oc n. In the case of elastic Maxwell molecules, vq is independent of temperature. The extension of this feature to 
the inelastic case defines Model A. On the other hand, one can assume that Vq is an increasing function of temperature 
(Model B). In particular vq(T) oc nT 1 / 2 makes Model B mimic the properties of IHS. 



III. UNIFORM SHEAR FLOW. RHEOLOGICAL PROPERTIES 



A. Hierarchy of moment equations 

Let us assume that the gas is under the USF. As said in the Introduction, this state is macroscopically defined by a 
constant density n, a spatially uniform temperature T(t), and a linear flow velocity u(y) = ayx [4j|. At a microscopic 
level, the USF is characterized by a velocity distribution function that becomes uniform in the local Lagrangian frame, 
i.e., 

/(r,v;t)=/(V,t). (34) 
In this frame, the Boltzmann equation (fTJ) reduces to 

d t f(V)-aV y ^rf(V) = J[V|/,/]. (35) 
Equation (|35|) is invariant under the transformations 

(V X ,Vy) -> (-V x ,-Vy), (36) 



Vj^-Vj, j?x,y. (37) 

This implies that if the initial state /(V, 0) is consistent with the symmetry properties (|36[) and (|3T|) so is the solution 
to Eq. (l35|) at any time t > 0. Even if one starts from an initial condition inconsistent with (|36|) and ((37|) . it is 
expected that the solution asymptotically tends for long times to a function compatible with (|36p and (|37|) . The 
investigation of this expectation, at the level of moments of degree less than or equal to four, is one of the objectives 
of this paper. 

The properties of uniform temperature and constant density and shear rate are enforced in computer simulations 
by applying the Lees-Edwards boundary conditions [13, HH , regardless of the particular interaction model considered. 
In the case of boundary conditions representing realistic plates in relative motion, the corresponding nonequilibrium 
state is the so-called Couette flow, where density, temperature, and shear rate are no longer uniform [54| . 

Multiplying both sides of Eq. by Y r u li2 _ is (y) and integrating over V one gets 

dtM-r\iii t 2...i s ~t~ ^ N r \i-^i 2 . . .i s J r\i\i2 ■ ■ -i s i (38) 

where we have called 

dV f(V)V y —Y rlili2 ... is (V). (39) 

In particular, 

2 1 

iY 2 |o = 2M \ xy , N \ yy — —-M \ xy , N \ xy = M \ y y + -M 2 \q, (40) 
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1 2 

No\ij = M \i y Sj x + M y y 5 ix + -M 2 \o(5i X 5jy + S JX 5 iy ) — -M Q \ xy 8ij, (41) 

A^4|o = 4M 2 | xy . (42) 

More in general, since V y dv,Y r \ ili2 is (y) is a polynomial of degree 2r + s, the quantity N r \ ili2 _ is can be expressed as 
a linear combination of moments of the same degree. In addition, thanks to the structure of the collision operator for 
IMM, the collisional moments J r u 1 , 2 ...i s only involve moments of degree equal to or smaller than 2r + s. Consequently, 
the hierarchy (|38p can be exactly solved in a recursive way. We will call asymmetric moments those that vanish for 
velocity distributions compatible with the invariance properties (|36l) and (|37p . The remaining moments will be referred 
to as symmetric moments. In particular, all the moments of odd degree are asymmetric. Among the moments of even 
degree, M 2r \ xz and M 2r \ xxxy , for instance, are also asymmetric. 

In the particular case of the moment M 2 \o = nTd/m, Eq. (|38[) becomes 

d t M 2]0 + 2aM 0lxy = -CM 2 | , (43) 

where use has been made of Eq. (JT7J) . This is not but the balance equation for the energy in the USF. It is convenient 
to introduce the scaled moments 



1 2T 2M 2l0 

M 2r|ili2 .., a = -^M 2r]ili2 .., s , „ = W_ = y_, (44) 



vq being the thermal speed, and the reduced shear rate 



a* = — . (45) 
vq 

In terms of these scaled variables, Eq. (f38|) can be rewritten as 

^^..i, + O^rVfa-i- " ( r + |) (<* + ^ a * M 0\ X y) K\i li2 ...i s = K\i li2 ...i 3 > (46) 

where £* = C/ v q and 

-^r|iii 2 ...i s = 2r+s^ r \W2— *»' "^rlf^...** = 2r+s ^r|ii i 2 ...i a ■ (47) 

It is apparent that the evolution equation (j46|) involves the second-degree moment ^olxy = ^xy/^P, which is the 
(reduced) shear stress. This quantity, along with the normal stress differences Mq, = (P xx —p)/2p and M^\yy = 
(P yy —p)/2p, are the most relevant ones from a rheological point of view. They will be analyzed in the next subsection. 



B. Second-degree moments 
From Eq. (j46|) one gets a coupled set of equations for the moments and M^ yy : 

^d t M* lxy + a* (M^ yy + + L 0{2 - 2«*MS\*y) M Q\*y = 0, (48) 



v 

where we have introduced the (reduced) shifted quantities 



dtM* lyy - -,a*M* Q]xy + u; | 2 - -,a*M*\ xy M* ]yy = 0, (49) 



so that 



"2r\s ~ {r + s/2)C 

w 2r \ s = , (50) 

vq 



W |2 = -. ■ (51) 



To close the set, we need in general the evolution equation for the reduced shear rate a*. From Eq. (|43p is straight- 
forward to obtain 

-d t a =a (C (52) 
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1. Model A. Hydrodynamic solution 



In Model A the collision frequency Vq is independent of temperature and thus it is a constant. Consequently, 
d t a* = so that the reduced shear rate a* remains in its initial value (regardless of the value of the coefficient 
of restitution a) and represents a control parameter measuring the departure of the system from the homogeneous 
cooling state. 

As in the elastic case [U it is easy to check that, after a certain kinetic regime lasting a few collision times, 
the scaled moments Mq\ xv and MZ reach well-defined stationary values, which are nonlinear functions of a and 



a* = a/vQ. From Eqs. 



and (|4*9")) . one has 



M, 



Q\xy 



1 - r M ^y 



= a \ M, 



l 0\yy 



(53) 



M o\ yy I 1 



^ a ^0\xy 



(54) 



for the stationary values, where we have defined 



4a* 



u i 2 (1 + a) 5 



(55) 



The solution to the set of equations (|53|) and (|54l) is 

7(a) 



0\yy 



M* 



dl{a) 



1 + 27(a) ' 



2 a 



a/2 



[l + 27(o)r 



(56) 



where 



7(a) = 3 sinh 



- cosh" 1 | 1 
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27, 



(57) 



is the real root of the cubic equation 



7 (l + 2 7 ) 2 = — . 



(58) 



Note that the reduced second-degree moments depend on a and a* through the scaled quantity a only. From Eq. 
([4"6"| it is also easy to prove that, for long times, the normal stresses M* z , . . . , M£ d along directions orthogonal to the 
shear plane xy are equal to M* y . Consequently, Mt x = —{d— l)MZ y . Analogously, the asymmetric second-degree 
moments (i.e., all the off-diagonal elements M^ - except M^ xy ) vanish. 

It is convenient to define a nonlinear shear viscosity rf and a viscometric function \I/ as 



P xy 



rf( a *) = _-2i£E = _ 2 

pa a 



0\xy 



(59) 



From Eqs. ([56")) and we have 



*(o*) 



v 2 p _ p M* - M*. 

t/Q r xx r yy r) \j\xx 0\yy 



p a* 



(60) 



r)*(a*) 



l + aj [1 + 2 7 (S)] 2 ' 



(61) 



*(a*) 



1 + aJ [l + 2 1 (a)Y 



(62) 
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Interestingly enough, the combination 



2 

a 



(63) 



is independent of the shear rate. Moreover, in the limit of small shear rate (for fixed a), Eq. (fS"8")) implies that 7^0, 
so that Eqs. (fBTj) and (f6"!?|) reduce to 

4 32 
{1+a.y (1 + a) 4 

The quantities -q* (0) and ^(O) are the NS shear viscosity and the Burnett value of the viscometric function, respectively, 
of Model A. 

It is important to remark that, although the scaled moments reach stationary values, the system is not in general 
in a steady state since the temperature changes in time. Actually, inserting the second expression of (|56p into Eq. 
(|43| . we get 

-a t lnr = -2w | 2 [7.-7(a)]. (65) 

where we have called 

C * d + 21 - a . (66) 



2w | 2 2d l + a 

Equation (|65|) shows that T(t) either grows or decays exponentially. The first situation occurs if 7(a) > 7 S - In that 
case, the imposed shear rate is sufficiently large (or the inelasticity is sufficiently low) as to make the viscous heating 
effect dominate over the inelastic cooling effect. The opposite happens if 7(a) < 7 S . A perfect balance between both 
effects takes place if 7(a) = j s . Inserting this condition into Eq. (|58|) one gets the steady-state point 



In this state, the a-dependence of the rheological properties is 

tf = »7*(0=(VrT ) ' ^ = *«) = T^f7TT V- ( 68 ) 

\a + 1 — a J 1 + a \a + 1 — a J 

Equations (|67p and (f6"8"]l agree with the results reported in Ref . [37| , while the more general expressions (fBTj) and ([6^)) 
had not been previously derived. 

Figure [1] shows the shear-rate dependence of the rheological functions f]*(a*) and ^(a*) for d = 3 and three values 
of the coefficient of restitution: a = 0.5 (highly inelastic system), a = 0.7 (moderately inelastic system), and a = 1 
(elastic system). The steady-state points (a*, 77*) and (a*, are also represented by circles for each one of the values 
of a. Given a value of a, the steady-state point splits each curve into two branches: the one with a* > a* corresponds 
to 7(a) > 7 S and so the temperature increases in time, while the branch with a* < a* corresponds to states with 
a decreasing temperature. We observe that, for a given value of a, both r\* (a*) and ^(a*) decrease with increasing 
shear rate. In the region of high shear rates (say a* > 1.5), the curves are practically insensitive to the value of the 
coefficient of restitution. For small and moderate shear rates, however, the influence of a is noticeable: at a given 
value of the reduced shear rate a*, the rheological quantities rf{a*) and \I/(a*) increase as a decreases. On the other 
hand, the steady-state values 77* and ^>* s decrease with increasing dissipation. 



2. Model B. Steady-state solution 



In Model B the collision frequency vq(T) is an increasing function of temperature, and so the reduced shear rate 
a* is not constant. The corresponding steady-state solution is obtained from Eqs. (|43|) - ([52|) by setting dt — > 0. It is 
given again by Eqs. (|56"l) - (f58")) . except that now 7(a) — > j s and a —* a*/w |2) where 7 S and a* are given by Eqs. (|6"6"]) 
and (|67[) . respectively. Therefore, the steady-state results are "universal" in the sense that they hold both for Model 
A and Model B, regardless of the precise dependence vq{T). 
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FIG. 1: (Color online) Plot of rj*(a*) (top panel) and ty(a*) (bottom panel) as functions of a* for d = 3 and a = 0.5 (solid 
lines), a = 0.7 (dashed lines), and a — 1 (dotted lines). The dash-dotted lines are the loci of steady-state points (a*,r/*) and 
(a*, ^s). They intercept the curves representing rj*(a*) and ty(a*) at the steady-state values indicated by circles. Note that 
the loci end at the points (aj,?7s) = (1.054, 0.563) and (aj, 9 a ) = (1.054, 1.688) corresponding to a = (represented by filled 
circles) . 



In order to have M*^ xy (a*) and M*^ yy (a*) in Model B, one has to solve numerically the nonlinear coupled set 
discard the kinetic stage of the evolution, and eliminate time in favor of a* [3]. The resulting rheological 
curves are illustrated in Fig. [2]at a — 0.7 and for the choice i>o{t) oc T 1 / 2 in Model B. Comparison with the analytical 
results corresponding to Model A shows that the influence of the temperature dependence of vq on the rheological 
properties is only significant for reduced shear rates smaller than the steady-state one. Since in this paper we want to 
focus on analytical results, henceforth we will only consider Model A, except for what concerns the common steady 
state of Models A and B. 



IV. THIRD- AND FOURTH-DEGREE MOMENTS 



In this Section we will analyze, in the context of Model A, the time evolution and the stationary values of the 
(scaled) third- and fourth-degree moments in the USF problem. The results will depend on both the reduced shear 
rate a* and the coefficient of restitution a, while they only depend on the latter in the common steady state. 

Let us assume that the scaled second-degree moments have reached their stationary values given by Eq. (|56[) . 
Therefore, Eq. (jit!]) becomes 

+ **K\i li2 ...i s -{ r + ^j [C - 2^o|27(S)] M r * |ili2 ... is = J r *| M2 ... ls , (69) 
where ds = v^Ai. In what follows we will particularize to a three-dimensional gas (d= 3). 
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FIG. 2: (Color online) Plot of rf(a*) (top panel) and ^(a*) (bottom panel) as functions of a* for d — 3 and a — 0.7. The 
dashed and solid lines correspond to Model A [Eqs. (|61|l and (|61|) ] and Model B with vo{T) oc T 1 ^ 2 [numerical solution of Eqs. 
(|48l) - (|52p ]. respectively. The circles represent the steady-state points (a*,rj*) and (a*, 9 S ), which are common to Models A and 



B. 



A. Third-degree moments 



As said in the preceding Section, all the third-degree moments are asymmetric and so they are expected to vanish 
for long times. Here we want to confirm this expectation and get the corresponding relaxation rates. 
In a three-dimensional system, there are 10 independent third-degree moments. Here we take 



V^2|x' ^2\y' M 2\z' M 0\xxy' M 0\xxz' M 0\xyy' M 0\yyz' ^0\xzz> M%*z> M 0\xyz\- 

From Eq. (|rj5|) and making use of the third-degree collisional moments (|18p and of the definition 
following set of equations: 



[9 S + ^o|3 + 3w | 2 7(a)] 



_ A/f * 
^ lvl 0\xxz 



M, 



0\yyz 



= o, 



(70) 

one gets the 
(71) 



[d s + u \3 + 3^ |27(a)] ( j M o\xx 



M, 



0\yzz 



0, 



(72) 



d s +u> 2 \i +3cj |27(a) 

d s 

i „* 





- ^o|3 + 3wo| 2 7(a) 
a* 



2a* 

2 



d s + lu \ 3 + 3cj | 2 7(a 




0\xyz 



(73) 
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FIG. 3: (Color online) Plot of the smallest eigenvalue, ^ m i n , associated with the time evolution of the third-degree moments as 
a function of a* for d = 3 and a = 0.5 (solid line), a = 0.7 (dashed line), and a = 1 (dotted line). 



la* 



I d s +uj 2 \i + 3w |27( a ) 

V t 

M* 

0\xxy 

V M* J 

\ 0\xyy / 



U* 2a* 

5 ^ 

d s + uj 2 \i + 3w |27( a ) 

<9 S + w |3 + 3w | 2 7(a) 



= o* I jM*,^ + M* lyzz 




20 " 





2a* 
fa* 

5 

d s + Wo|3 + 3cj | 2 7(a) / 



(74) 



[<9 S + w | 3 + 3w | 27 (a)] M*|^ z = a* 



2 2 

5 M 2*| a + 5 M 0%a;y _ M 0\yzz 



The characteristic equations associated with Eqs. ATI"]) - ([75]) are 



L = ^o|3 + 3cj | 2 7(a) = -{I + a)' 



7(a) + 2 



(75) 



(76) 



[cj | 3 + 3w | 2 7(a) - (] 2 [w 2 |i + 3cj | 2 7(a) - £] = 



2(w | E 



" w 2|lJ *2 

—a , 



(77) 



[w |3 + 3w |27( a ) - [ w 2|i + 3w | 2 7(a) - £]' 



2(a; | 3 - a; 2 |i] 
25 



a* 2 [7^ |3 + 23^,1 + 30 (3^ | 2 7(a) - l)] , (78) 



where £ denotes the corresponding eigenvalues. The time evolution for long times is governed by the eigenvalue £ m in 
with the smallest real part. It can be checked that £ m [ n [which is the smallest real root of the quartic equation (|78p ] 
is positive definite for all a and a*. Consequently, all the scaled third-degree moments vanish in the long time limit, 
as expected by symmetry arguments. In addition, at a given value of a*, the larger the inelasticity the longer the 
relaxation time (which is proportional to ^~j n ). The shear-rate dependence of £ m - ln is plotted in Fig. [3] for the same 
values of a as considered before. It is interesting to remark that £ m j n is not a monotonic function of a*, reaching a 
minimum value at a certain shear rate. 



B. Fourth-degree moments 



In a three-dimensional system, there are 15 independent fourth-degree moments, of which 9 are symmetric and 6 
are symmetric, in the sense described at the beginning of Sec. IIIII 
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1. Asymmetric moments 



Because of the symmetries of Eq. (|35f . the symmetric and asymmetric moments are uncoupled. Although the 
relevant moments are the symmetric ones, we first analyze the time evolution of the asymmetric moments, for the 
sake of completeness. As the set of asymmetric moments, we choose 



i-^2\xz' M 2 \ yx , M \ xxxzy , M Q \ yyyz , M Q \ XZZZ , M \ yzzz }. 
They verify the following set of equations: 

[d s + cj |4 + 4cj | 2 7(a)] {M^ yyyz ~ M* ]yzzz ) = 0, 



(79) 
(80) 



/ d s + lo 2 \ 2 + 4w | 2 7(a) 





49 



d s +uj 2 \2 + 4w |27(a) 





-2a* 



-2a* 




-±§a* d s + w |4 + 4w | 27 (a) 

fa* d s + uj 014 + 4^ | 27 (S) / 



-Ala* 

14" 



M. 



M. 



2\j 



!)!)!)■- 



M, 



M 0\xxxz + M 0\xzzz 



-a* (MS 



— M* ) 

\yyyz lv± 0\yzzz) 



(81) 



Q\yzzz 



[d s + ^ |4 + 4^ | 27 (a)] (M* lxxxz - M* lxzzz ) = -a*(M* lyyyz + M* {yzzz ) - -a*(M* lyyyz - M* lyzzz ). (82) 
The eigenvalues associated with the time behavior of the asymmetric fourth-degree moments are 

I = w | 4 + 4w | 2 7(a) = (1 + a) 2 
and the roots of the characteristic quartic equation 

[won + 4 ^o|27(S) " e] 2 [W2|2 + 4^ |27(«) - t\* = ^"'^g — ^ [3^ | 4 + llw 2 | 2 + 14 (4^ | 4 7(a) - l)] ■ (84) 

All the eigenvalues have a positive real part. Therefore, all the asymmetric moments defined in Eq. (|79[) decay to zero 
in the long time limit, the final stage being characterized by the smallest real root £ m [ n of Eq. ([84]) . This eigenvalue 
is plotted in Fig. |4] for the same cases as in Fig. [3] Again, a non-monotonic behavior is observed. On the other hand, 
for given values of a* and a, the value of £ min corresponding to the third-degree moments is smaller than the one 
corresponding to the asymmetric fourth-degree moments. 



7(«) 



50 + 7a 



126 



(83) 
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2. Symmetric moments 
In parallel to the elastic case [47|, we choose the following set of 9 symmetric moments 

V^4|0> ^2 1 xx' M 2\yy> M 2\xy> M 0\xxxxi M 0\yyyy M 0\zzzz' M 



M*, \ 

xxxy 0\xyyy S ' 



The combination 



M 9 = 3MZ 



Olxxxx 



AM, 



0\yyyy 



4K\zzzz = "4 / dV (6V*V Z 2 K K 4 ) /(V) 



is the average of a quantity independent of V x , so the associated combination 3N^ XXXX — 47V,* 
Moreover, the combination 3 JZ — 4J*, — 4 J*, 

' 0\xxxx 0\yyyy Q\zzzz 

Eq. (H yields 



[<9 S + cj | 4 + 4w | 27 (a)] X 9 = 0. 

The solution to this equation is simply 

Mg(a) = M 9 (0)e- e ° s , £ g = u 0l4 + 4^ | 27 (u). 

Since cj |4 > 44], one has £g > and so Mg decays to 0. 

The remaining eight moments in (|85p are coupled. In matrix form, Eq. (|69p becomes 



(85) 

(86) 
vanishes. 

j/ |4jVfg due to the fact that M ^ yy = M \ zz . Therefore, 

(87) 



\yyyy 



■4N* 



{8^,d s + C aal )M c 
Here, M. is a vector made of the following 8 moments 



Ccr, cr=l,. 



M = 



( M* m \ 

M. 
M. 
M, 



2\xx 

■* 

2\yy 



0\yyyy 
1V1 0\zzzz 



M. 
M* 



2\xy 



and the square matrix C is 



where X is the 8x8 identity matrix and 



Q\xxxy 
V M 0\xyyy J 



C = 4w i 2 7(a)^ + 



C' = 
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<^4|0 



















^2|2 



















<x>2|2 



















^0|4 



















w 0|4 




^a* 




^a* 








15 " 


fa* 




-|a* 


-|a* 







49 " 


49" 


- 5 a* 


-14«* 


V 





6_ * 

49 " 


i^a* 
49 u 


2a* 


ia* 
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\ 


32 * 
2l" 

-Pa* 
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2a* 




_- ii a * 

^a* 
245" 




fa* 


-fa* 
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^2|2 













^0|4 













w 0|4 


/ 



In addition, C is a vector of elements made of second-degree moments, namely 

Ci - \K - K (m* 0V J + M* [yy 2 + M* { J + 2M^ xy 2 ) , 



-A4 (^A/,*^ 2 



M, 



0\yy 



M, 



0\zz 



M, 



0\xy 



(88) 



(89) 



(90) 



(91) 



(92) 



(93) 
(94) 
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FIG. 5: (Color online) Plot of the smallest eigenvalue, l-mi-n, associated with the time evolution of the symmetric fourth-degree 
moments as a function of a* for a — 0.5 (solid line), a — 0.7 (dashed line), and a = 1 (dotted line). The circles indicate the 
location of the corresponding values of the critical shear rate. 



C 3 = 2 A 3 M o, ra 



M, 



51M, 



6M 



0|:rx 



Ola; a; 



OI2Z 



M, 



0|a;y 



6M n *, 



0|zz 



48M, 



0|a;y 



(95) 
(96) 



C 5 = ^rAj ( 51M, 



6Mn 



6M, 



12M, 



(97) 



C R = -ASM, 



0|a:y 



(98) 



c 7 - =K M o\x 



2M, 



0\yy 



(99) 



c H = =a;m 



0|a;y 



5M 



2M 



The solution of Eq. 



where 



can be written as 
M(s) 



-Cs 



[M(0) - Moc] + M c 



cr x c. 



(100) 



(101) 



(102) 



Similarly to the cases discussed above, the long time behavior of M.„ (a = 1, . . . , 8) is governed by the eigenvalue 
^min of the matrix C with the smallest real part. We have checked that £ m in is a real quantity that, for a given value 
of a, monotonically decreases with increasing shear rate. It is plotted in Fig. O for a — 0.5, 0.7, and 1. The most 
important feature of Fig. [5] is that, for any given value of a, £ m i n becomes negative for shear rates larger than a certain 
"critical" value a* (a). This means that, if a* > a*, the symmetric fourth-degree moments exponentially grow in time. 
This singular behavior of the scaled moments implies that the velocity distribution function (scaled with the thermal 
speed) develops an algebraic high-velocity tail in the long time limit. It is interesting to remark that this effect is 
also present in the elastic limit, where it has been extensively studied [13, IH, |49[. As observed in Fig. [5l the main 
influence of inelasticity is to decrease the value of the critical shear rate. 

Let us analyze the phase diagram associated with the singular behavior of the fourth-degree moments. This is 
shown in Fig. O where the curve a* (a) splits the parameter space into two regions: the region below the curve 
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a 



FIG. 6: (Color online) Phase diagram for the asymptotic long time behavior of the fourth-degree (symmetric) moments. The 
shaded region below the curve a? (a) (thick solid line) corresponds to states with finite asymptotic values of the scaled fourth- 
degree moments, while the region above the curve defines the states where those moments diverge in time. The dash-dotted 
line represents the steady-state points al(a). It intercepts the critical curve at(a) at the point (a, a*) — (0.046, 1.041). 




1 LJ i 1 ■ 1 ' lJ ' i 1 
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a 

FIG. 7: (Color online) Plot of the asymptotic long time value of the scaled moment M^ as a function of a* for a = 0.5 (solid 
line), a — 0.7 (dashed line), and a — 1 (dotted line). The dash-dotted line represents the values of M^ at the steady states 
aj(a) for 0.046 < a < 1. It intercepts the curves representing M^ (a*) at the points indicated by circles. The vertical dotted 
lines are the asymptotes of the curves. 

corresponds to states (a, a*) with finite asymptotic values of the scaled fourth-degree moments (i.e., £ m - ln > 0), while 
the region above the curve defines the states where those moments diverge in time {t m \ n < 0). Figure [5] also includes 
the locus of steady-state points (a, a*) [cf. Eq. (|67p ]. Below the latter curve the inelastic cooling dominates over 
the viscous heating and so the temperature decreases in time, while the opposite happens above it. It is apparent 
that the curve a*(a) lies inside the region a* < a*(a), except for the small interval a < 0.046 or, equivalently, 
a* (0.046) = 1.041 < a* < a*(0) = 1.054. In conclusion, in order to find diverging moments, one has to consider states 
with rather large values of the shear rate at which the viscous heating is much higher than the collisional cooling. 

For states with a* < a* (a) the scaled (symmetric) fourth-degree moments reach well-defined finite values in the 
asymptotic long time limit. From Eq. (|101[) one has 

lim M(s) = Mao, (103) 

s — >oo 

where Aioo is defined by Eq. (|102p . As an illustration, Fig. [7] shows the shear-rate dependence of the asymptotic long 
time values of M* = (n^)" 1 / dW 4 f(V) for a = 0.5, 0.7, and 1. The values of M* |0 at a* = correspond to the 
homogeneous cooling state [3]. We observe that, given a value of a, the scaled moment M^ rapidly increases with 
the shear rate, having a vertical asymptote at a* — a* (a). Moreover, for a given value of a* , the value of the moment 
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increases with dissipation. Figure [7] also includes the curve representing the values of M£, Q at the steady states a*(a). 
This curve has a vertical asymptote at a* — 1.041. 

V. CONCLUDING REMARKS 

The simple or uniform shear flow is perhaps the most widely studied inhomogcneous state for elastic and inelastic 
gases. Despite its apparent simplicity, this state has proven to be useful to shed light on the nonlinear response 
of the system to the presence of strong shearing. This response is conventionally measured by the non-Newtonian 
rhcological properties (namely, the nonlinear shear viscosity r\ and viscomctric function ty), which arc related to the 
second-degree velocity moments (pressure tensor). On the other hand, higher degree moments are also important to 
provide indirect information on the features of the velocity distribution function. For inelastic gases, there are two 
relevant control parameters: the shear rate a scaled with an effective collision frequency vq, i.e., the reduced shear 
rate a* = a/vg, and the coefficient of normal restitution a. When the velocity moments are conveniently scaled with 
the thermal speed, they are expected to become, after a kinetic transient regime, nonlinear functions of both control 
parameters. 

To address the above issues in the context of the Boltzmann equation without having to resort to approximate 
methods or computer simulations, one can consider simplified collision models. As in the case of elastic collisions 
[47l | . the inelastic Maxwell model (IMM) renders itself to an analytical treatment. Here, taking advantage of a recent 
derivation by the authors of the collisional moments for IMM through fourth degree |44|, we have determined the 
pressure tensor and the fourth-degree velocity moments of the USF problem in an exact way. Two different classes of 
IMM have been defined: Model A, where v is independent of temperature, and Model B, where v$ is an increasing 
function of temperature. In the USF the temperature changes in time due to two opposite effects: viscous heating 
and collisional cooling. Therefore, a steady state is eventually achieved in Model B when both effects cancel each 
other. However, in Model A a steady state does not generally exist, except for a specific value a* = a* s (a), Eq. (j6"Tj) . 
It is important to note that the results in the steady state are the same for both classes of models. Since the reduced 
shear rate a* changes with time in Model B, the goal of obtaining the velocity moments as functions of a* and a 
requires the use of numerical tools in that case. However, a* — const in Model A and so the independent influence of 
a* and a can be studied analytically. Thus, in this paper we have focused on Model A, except in what concerns the 
steady state which, as said before, is common to Models A and B. 

As mentioned above, the relevant transport properties in the USF problem are the nonlinear shear viscosity rj* (a*) 
and the viscometric function '5 (a*). Their explicit forms are given by Eqs. (|B"Tj) and (|62|) . respectively. These results 
extend to inelastic collisions the expressions obtained long time ago by Ikenberry and Truesdell for (elastic) Maxwell 
molecules [51(. With respect to the dependence of rj*(a*) and \l/(a*) on inelasticity, our results show that its influence 
is quite significant for small and moderate shear rates, both rheological properties being practically insensitive to 
dissipation in the region of high shear rates. In the steady-state solution, rj*(a*) and \I/ S (a*) decrease when decreasing 
the coefficient of restitution. Moreover, as expected, the (scaled) third- and asymmetric fourth-degree moments 
vanish in the long time limit. Consequently, beyond the rheological properties, the next nontrivial moments are the 
symmetric fourth-degree moments. An important result is that, for a given value of the coefficient of restitution, 
these moments are divergent for shear rates larger than a certain critical value a* (a). This singular behavior is also 
present in elastic systems [13, EH [4!| , where it has been shown that this divergence is consistent with an algebraic 
high- velocity tail of the velocity distribution function. The main effect of inelasticity is to decrease the value of a* (a) 
as the gas becomes more inelastic. In addition, the phase diagram associated with this singular behavior shows that 
the value of a* (a) is rather large in the whole domain < a < 1, so that in order to get diverging moments one has to 
consider states at which the collisional cooling is strongly dominated by the viscous heating effect. As a consequence, 
nonlinear shearing effects are still significant for a* < a*, as illustrated in Fig. [7] for the scaled moment M*^ Q . 

The results derived in this paper can be useful for analyzing different situations. First, the knowledge of the 
shear-rate dependence of the second- and fourth-degree moments of the USF allows one to determine the generalized 
transport coefficients characterizing transport around the simple shear flow [461 ] . Another possible direction of study 
is the extension of the present analysis for the rheological properties to multicomponent systems. Previous works 
carried out for IMM [3?], [39[ have shown the tractability of the Maxwell kinetic theory for these complex systems and 
stimulate the performance of this study in the near future. 
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